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^ ; ABSTRACT 

o 

Observations of dark matter dominated dwarf and low surface brightness disk galaxies favor density 
£S) \ profiles with a flat-density core, while cold dark matter (CDM) N-body simulations form halos with 

central cusps, instead. This apparent discrepancy has motivated a re-examination of the microscopic 
' nature of the dark matter in order to explain the observed halo profiles, including the suggestion that 

CDM has a non-gravitational self-interaction. We study the formation and evolution of self-interacting 
dark matter (SIDM) halos. We find analytical, fully cosmological similarity solutions for their dynamics, 
which take proper account of the collisional interaction of SIDM particles, based on a fluid approximation 
derived from the Boltzmann equation. The SIDM particles scatter each other elastically, which results in 
an effective thermal conductivity that heats the halo core and flattens its density profile. These similarity 
solutions are relevant to galactic and cluster halo formation in the CDM model. We assume that the local 
density maximum which serves as the progenitor of the halo has an initial mass profile SM/M oc M~ £ , 
as in the familiar secondary infall model. If e = 1/6, SIDM halos will evolve self-similarly, with a cold, 
supersonic infall which is terminated by a strong accretion shock. Different solutions arise for different 
values of the dimensionless collisionality parameter, Q = o~pbr s , where a is the SIDM particle scattering 
cross section per unit mass, pb is the cosmic mean density, and r s is the shock radius. For all these 
| solutions, a flat-density, isothermal core is present which grows in size as a fixed fraction of r s . We find 

. two different regimes for these solutions: 1) for Q < Qth{— 7.35 x 10~ 4 ), the core density decreases and 

core size increases as Q increases; 2) for Q > Qth, the core density increases and core size decreases as 
Q increases. Our similarity solutions are in good agreement with previous results of N-body simulation 
of SIDM halos, which correspond to the \ow-Q regime, for which SIDM halo profiles match the observed 
galactic rotation curves if Q ~ [8.4 x 10~ 4 — 4.9 x 10~ 2 }Qth, or a ~ [0.56 — 5.6] cm 2 g~ 1 . These similarity 
solutions also show that, as Q — > oo, the central density acquires a singular profile, in agreement with 
some earlier simulation results which approximated the effects of SIDM collisionality by considering an 
ordinary fluid without conductivity, i.e. the limit of mean free path \ m f p — > 0. The intermediate regime 
where Q ~ [18.6 — 231}Q t h or a ~ [1.2 x 10 4 — 2.7 x 10 4 ] cm 2 g~ 1 , for which we find flat-density cores 
comparable to those of the low-Q solutions preferred to make SIDM halos match halo observations, has 
not previously been identified. Further study of this regime is warranted. 

b : 
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I. INTRODUCTION 

The cold dark matter (CDM) model provides a suc- 
cessful framework for understanding the formation and 
evolution of structure in the universe. According to this 
model, structure forms out of primordial density fluc- 
tuations by hierarchical clustering; small objects form 
first, and then merge to make bigger objects. However, 
this CDM model has several problems. 

There has been a recent concern about the possible 
discrepancy in the inner density profile of dark-matter 
dominated cosmological halos between N-body simula- 
tions and observations. CDM halos in N-body simula- 
tions have a density cusp (p oc where (3 ranges from 
-1 (Navarro, Frenk & White 1997; hereafter "NFW") 
to -1.5 (Moore et al. 1999; hereafter "Moore profile"), 



while observations of dwarf and LSB galaxies, which 
are believed to be dark-matter dominated, indicate flat- 
density (soft) cores. 

Spergel & Steinhardt (2000) suggested that the purely 
collisionless nature of CDM be replaced by self-interacting 
dark matter (SIDM) to resolve this discrepancy. The 
non-gravitational, microscopic interaction (e.g. elastic 
scattering) of SIDM particles is postulated to be strong 
enough to produce a soft core by heat conduction. Cos- 
mological N-body simulations which incorporate a fi- 
nite scattering cross-section a for the SIDM particles 
show that this scheme successfully produces soft cores 
in halos (e.g. Dave et al. 2001; Yoshida et al. 2000b). 

Previous analytical study (Balberg, Shapiro & In- 
agaki 2002) and estimates based on N-body simula- 
tions (Burkert 2000; Kochanek & White 2000) indicate 
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that SIDM cores will become unstable by undergoing 
gravothermal catastrophe in a Hubble time. This po- 
tential instability of SIDM cores means that, if the mat- 
ter content of the universe is really SIDM particles, we 
are living in a special epoch where SIDM halos still have 
soft cores. However, these analyses are based on iso- 
lated halos. In the context of the hierarchical clustering 
scenario, cosmological infall is inevitable. We are led 
to pose the following question: Can cosmological infall 
delay the core collapse of SIDM halos? 

In this paper we study the formation and evolution 
of SIDM halos with a proper treatment of cosmologi- 
cal infall. We find that, for a range of infall rates, the 
collapse (i.e. gravothermal catastrophe) of SIDM cores 
can be completely inhibited. We also find a new range 
of values of the scattering cross section which can pro- 
duce soft cores, separate from and in addition to the 
range previously identified by N-body simulations. To- 
ward this end, we find a self-similar solution to describe 
the evolution of SIDM halos for arbitrary degree of col- 
lisionality which allows a detailed analytic study. 

II. BASIC EQUATIONS 

We use a fluid approximation derived from the Boltz- 
mann equation to describe the dynamics of SIDM halos. 
The basic assumptions we make are as follows: 1) the 
system has spherical symmetry; 2) infall is radial, cold, 
and continuous; 3) virialized objects have an isotropic 
velocity distribution. The first two assumptions have 
been used before for secondary infall models (e.g. Fill- 
more & Goldreich 1984). The third assumption comes 
from the property of CDM halos in N-body simulations: 
the velocity distribution is isotropic at the center (i.e. 
radial velocity dispersion is equal to tangential velocity 
dispersion) and acquires only moderate anisotropy as 
one goes out to the virial radius. 

With these assumptions, the moments of the col- 
lisional Boltzmann equation (i.e. valid for arbitrary 
degree of collisionality) can be shown to reduce to the 
usual fluid equations for an ideal gas with 7 = 5/3, 
where 7 is the ratio of specific heats. This approach 
has been widely used in the literature of stellar dy- 
namics (e.g. Bettwieser 1983), especially for the study 
of gravothermal catastrophe. We note that the par- 
ticle components in both cases, that of stellar systems 
(stars) and of SIDM halos (SIDM particles), are weakly 
collisional, and a corresponding conductive heating 
term enters the energy equation in both cases, but 
it has a different origin in each case. For the stel- 
lar systems (number of particles ~ 10 4 — 10 6 ), weak, 
gravitational scattering is responsible for the collisions, 
while for SIDM halos, nongravitational, microscopic 
self-interaction is responsible for the collisions. 

In summary, we use the following fluid-like equa- 
tions of conservation of mass, momentum, and energy, 
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where p, u{= (v r )), p(= p{v r — u) ), and / are mass 
density, bulk radial velocity, pressure, and heat flux, 
respectively, M is the integrated mass inside radius r, 
and -§t = ■§[ + u ~§f- For a collisionless system, / = 
and we will call it an "adiabatic" process. For SIDM 
halos, heat conduction takes the following form: 
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where a and b are constants of order of unity (Balberg, 
Shapiro & Inagaki 2002). We consider an elastic scat- 
tering case, in which a = 2.26 and b = 1.002. 

III. SELF-SIMILAR EVOLUTION OF SIDM 
HALOS 

We seek a self-similar solution to describe the evo- 
lution of SIDM halos in an Einstein-de Sitter universe. 
This approach makes it possible to analyze properties 
of the solution in detail, because an exact solution can 
be found. Previous authors have found self-similar so- 
lutions for the evolution of cosmological halos for cold, 
collisionless dark matter in radial motion (e.g. Fill- 
more & Goldreich 1984) and for an ordinary, adiabatic, 
baryonic fluid (e.g. Bcrtschinger 1985). Our problem 
differs from these in that our gas is neither collision- 
less nor in purely radial motion and it differs from an 
ordinary, adiabatic, baryonic fluid, since our fluid-like 
conservation equations also contain the SIDM conduc- 
tion term. On the other hand, we were encouraged to 
look for a self-similar solution, since the effect of con- 
ductivity is to drive a gas towards isothermality, and 
there are self-similar equilibrium solutions known in the 
limit of isothermality which have soft cores, such as the 
"truncated isothermal sphere" (TIS) model proposed 
by Shapiro, Ilicv & Raga (1999). 

We start from an adiabatic (i.e. no conduction, or 
/ = 0), self-similar solution to find the condition re- 
quired for the self-similarity of SIDM halo solutions. If 
the initial ovcrdensity has a scale-free power-law form, 
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where e is any positive constant, then adiabatic infall 
will form a self-similarly growing object (Fillmore & 
Goldreich 1984). A cold, supersonic, radial infall leads 
to a strong, spherical shock which surrounds a subsonic 
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region of shock-heated gas. The corresponding length 
scale, e.g. a turnaround radius r ta , grows as 



r ta oc 



where 
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In this case, the shock always occurs at the same fixed 
fraction of r ta , the mean density inside the postshock 
region is a constant multiple of the time- varying cosmic 
mean density pb, and the radial profiles of all fluid vari- 
ables are self-similar. When heat conduction by SIDM 
particles is introduced, however, an additional length 
(time) scale is present which breaks the self-similarity 
of these adiabatic solutions, in general. 

There is one particular infall rate, however, which 
does not break the self-similarity even when heat con- 
duction is introduced. According to the adiabatic solu- 
tions, the thermal energy changes according to oc 
r 2 a t~ 5 , while the conductive heating term would intro- 
duce a change as 



d 

2 dr 



2 /) oc rf a t~ 7 . To preserve sim- 
ilarity after the addition of heat conduction, therefore, 
we must have r ta oc t 2 , or 
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Let us first examine properties of the adiabatic so- 
lution (/ = 0) for the case of e = 1/6. We especially 
focus on the density profile. The density profile of the 
adiabatic solution within the shock radius r s agrees 
quite well with the N-body results for CDM halos (see 
Fig. 1). The adiabatic solution is well-fit, for exam- 
ple, by an NFW profile with concentration parameter 

3 < cnfw < 4 (i.e. fractional deviation A < 17% for 
0.014 < r/r 2 oo < ^s/^200, where r 2 oo is the NFW pro- 
file radius within which (p) = 200pb, and r s ~ 0.6r2oo)- 
The inner density slope of the similarity solution for 

4 x 10~ 3 < r/r 20 o < 1-4 x 10~ 2 is -1.27, which is be- 
tween -1 (NFW) and -1.5 (Moore profile). A natural 
question is then: why does this adiabatic solution for 
e = 1/6 match the CDM N-body halos so well? 

The theory of structure formation from density peaks 
in the Gaussian random noise distribution of initial 
density fluctuations gives an interesting clue to this cor- 
respondence. According to Hoffman & Shaham (1985), 
local maxima of Gaussian random fluctuations in the 
density can serve as the progenitors of cosmological 
structures. For a power-law power spectrum of initial 
fluctuations, P(k) oc k n , the initial density profile can 
be approximated as y oc r K , where k = 3e = n + 3. 
The value e = 1/6 corresponds to k = 1/2, or n = —2.5, 
while n = —2.5 roughly corresponds to galaxy-mass 
structures in CDM, as shown in Fig. 2. 




Radius (i7r 200 ) 



Fig. 1. — Halo density profile for adiabatic solution with 
e = 1/6 for standard CDM halos, compared to NFW and 
Moore profiles. The adiabatic solution, whose inner slope 
is -1.3, is a good fit to density profiles of CDM halos from 
N-body simulations. 
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Fig. 2. — Effective index of the power spectrum of initial 
density fluctuations (P(k) oc k n "ff) vs. mass of halos at 
their typical formation epoch, for flat, cluster-normalized, 
matter-dominated CDM (SCDM) and COBE-normalized 
ACDM universes. 
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Fig. 3. — Dimensionless profiles for similarity solutions in the low-Q regime, where rt a is the time-varying turnaround 
radius and the shock radius is r s = 0.09rt a - Isothermal, flat-density cores exist. As Q increases, core density decreases and 
core temperature increases. Dimensionless similarity variables follow the definitions in Bertschinger (1985). 



V. SELF-SIMILAR SIDM HALOS (e = 1/6; 
WITH CONDUCTION) 

We find that different similarity solutions of SIDM 
halos arise for different values of the collisionality pa- 
rameter, Q = apbr s . The dimensionless Q is closely 
related to P = apv re iAt, the number of collisions each 
particle experiences during a time At in a local den- 
sity p with relative velocity v re i, sometimes used else- 
where to parameterize the collisionality of SIDM par- 
ticles (Dave et al. 2001). Roughly speaking, high Q 
means high collision rate. We also find that there are 
two different regimes: the low-Q regime and the high-Q 
regime. 

In the low-Q regime, the collision mean free path in 
the core is larger than the size of the halo. In this 
regime, the flattening of the core density profile in- 
creases as Q increases (see Fig. 3). As Q increases, 
that is, the central density is lower and the core radius 



is larger. 

In the high-Q regime, the mean free path in the core 
region is smaller than the size of the halo. The flatten- 
ing of the core density profile decreases as Q increases, 
because the mean free path decreases (diffusion limit). 
As Q increases in this regime, therefore, the core den- 
sity increases and the core radius shrinks (see Fig. 4). 
The limiting case of Q = oo corresponds to the solution 
with no conduction, which agrees with N-body simula- 
tions with infinite cross-section (Yoshida et al. 2000a; 
Moore et al. 2000). In this limit of maximal collisional- 
ity, the nonadiabatic solution approaches the adiabatic 
solution (see Figs. 3 and 4). 

An important aspect of this self-similar solution is 
that halo cores never enter the gravothermal catastro- 
phe phase, because cosmological infall causes kinetic 
energy to be pumped continuously into the core to pre- 
vent it. As a result, cores grow as a fixed fraction of the 
shock radius. As mentioned in Section IV, the particu- 
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Fig. 4. — Dimensionless profiles for similarity solutions in the high-Q regime. Profiles are indistinguishable from those in 
Fig. 3, even though Q values are quite different. The effect of Q is reversed from that of low-Q solutions: as Q increases, 
core density increases and core temperature decreases. 



lar value of e = 1/6 adopted here to make the evolution 
of SIDM halos self-similar (which, therefore, precludes 
core collapse), is natural for galaxy formation in a CDM 
universe, as can be understood in the context of the 
growth of density peaks. According to this theory, the 
average initial density profile responsible for the forma- 
tion of individual galaxy-mass halos in the CDM model 
corresponds to a value of n e ff which is not far from 
—2.5, the value required to make e = 1/6. As such, the 
self-similar behavior reported here for SIDM halos may 
be a realistic approximation for galactic halos. If the 
infall rate which builds an individual galactic halo even- 
tually tapers off at late times, however, as reported for 
N-body simulations by Wechsler et al. (2002), this will 
break the condition necessary for SIDM self-similarity 
and may allow gravothermal catastrophe to proceed. 



VI. THE ALLOWED RANGE OF SCATTER- 
ING CROSS - SECTION FOR AN SIDM 
UNIVERSE 

In a particular universe, we can relate a to Q if we 
know the typical formation epochs for halos of differ- 
ent masses, according to the Press-Schechter formalism. 
We find that a relatively narrow range of Q values, 

O.^-^x^lA)^), (9) 

characterizes the entire range of halo masses from 
dwarf-galaxy-mass to cluster-mass, in the currently- 
favored ACDM universe. We also find 

O , [ U-C,]x 1 0-»(^)(^ FT ) ,10) 

for a flat, cluster- normalized, matter-dominated CDM 
(i.e. SCDM) universe (See Fig. 5). 
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Fig. 5. — The collisionality parameter Q vs. mass of halos 
at their typical formation epoch for a = 5.6 cm 2 g^ 1 and 
h = 0.7. 

The range a — [0.56 — 5.6] cm 2 g _1 is the preferred 
range of the scattering cross section found by cosmo- 
logical N-body simulations for the ACDM universe to 
match observed galactic rotation curves (e.g. Dave et 
al. 2001). The above relation (equ. [9]) with h = 0.7 
then yields Q = [6.2 x 10~ 7 - 3.6 x 10" 5 ] for the ACDM 
universe, which is in the low-Q regime. The same pro- 
file shapes will occur in SCDM for the same Q- values. 
If we apply this low-Q solution to SCDM with h = 0.70, 
using equation (10), we get a = [0.23 — 3.3] cm 2 g~ 1 . 

As shown in Section V, there also exist high-Q solu- 
tions which yield profiles which are quite similar to the 
low-Q solutions which produce observationally accept- 
able soft cores. We find that Q = [1.37 x 10" 2 - 1.7 x 
10~ 4 ] in the high-Q regime produces profiles like those 
in the low-Q regime for Q = [6.2 x 10~ 7 - 3.6 x 10" 5 ]. 
From the relationship between a and Q (equs. [9] and 
[10]), we predict that a = [1.2 x 10 4 - 2.7 x 10 4 ] cm 2 ^ 1 
can also produce acceptable soft cores, therefore, for 
ACDM, while the corresponding a for SCDM will be 
about [5.11 x 10 3 - 1.56 x 10 4 ] cm 2 g-\ 

VII. CONCLUSION 

We have found similarity solutions which allow us to 
describe the dynamical origin and evolution of SIDM 
halos in a fully cosmological context, analytically, for 
the first time. Our solutions are based upon fluid-like 
conservation equations which can be derived by taking 
moments of the Boltzmann equation, assuming that 
the particle velocity distribution inside virialized ha- 



los is approximately isotropic, as CDM N-body simu- 
lations suggest. Our solutions confirm and explain the 
results of N-body simulations which have attempted 
to incorporate the collisional effects of SIDM numeri- 
cally within the CDM model and allow us to extend 
and generalize those results to a much wider range of 
parameters and scales than has so far been simulated. 
Along the way, we have also found a similarity solution 
for the adiabatic case of standard CDM (i.e. with no 
SIDM self-interaction) which serves to derive analyti- 
cally the halo profile of CDM halos found previously 
by N-body simulation (i.e. intermediate between NFW 
and Moore profiles), as a product of self-similar cosmo- 
logical infall. 

As seen in our similarity solutions, cosmological in- 
fall can affect the dynamics of virialization significantly. 
In the context of hierarchical clustering, in fact, the 
collapse of halo cores previously predicted for isolated 
SIDM halos is entirely prevented by cosmological in- 
fall, according to these similarity solutions. For realis- 
tic mass accretion histories in a CDM universe, there- 
fore, SIDM core collapse will be delayed until this in- 
fall becomes negligible. Accordingly, previous analy- 
ses (Burkert 2000; Balberg, Shapiro & Inagaki 2002) 
which predict core collapse in a Hubble time should be 
re-examined. 

By considering the full range of collisionality from 
Q = to Q = oo, we have discovered that the low-Q 
regime which has been found previously by N-body sim- 
ulation to produce observationally acceptable soft cores 
is not unique. We predict that a cross-section of the 
order of 10 4 cm 2 g~ x can also produce acceptable soft 
cores, in addition to the regime of smaller cross-section 
previously identified. Therefore, we suggest that cos- 
mological N-body simulations be performed which in- 
corporate SIDM particles with a = [5 x 10 3 — 5 x 
10 4 ] crr^g^ 1 to explore this possibility further. 

In reality, the mass accretion rates of CDM halos in 
N-body simulations are found to depart from the self- 
similar values over time, declining at late times rather 
than growing (Wechsler et al. 2002). We have accord- 
ingly begun to study the evolution of SIDM halos with 
a more realistic infall, using a 1-D, spherical hydrody- 
namics code, and will describe those results elsewhere 
(Ahn and Shapiro, in preparation). 

This work was partially supported by NASA ATP 
Grant NAG5-10825 and Texas Advanced Research Pro- 
gram Grant 3658-0624-1999. 
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